Magnetic properties of (VO)2P207: two-plane structure and spin-phonon interactions 
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Detailed experiments on single-crystal (VO)2P207 continue to reveal new and unexpected features. 
We show that a model composed of two, independent planes of spin chains with frustrated magnetic 
coupling is consistent with nuclear magnetic resonance and inelastic neutron scattering measure- 
ments. The pivotal role of PO4 groups in mediating intrachain exchange interactions explains both 
the presence of two chain types and their extreme sensitivity to certain lattice vibrations, which re- 
sults in the strong magnetoelastic coupling observed by light scattering. We compute the respective 
modifications of the spin and phonon dynamics due to this coupling, and illustrate their observable 
consequences on the phonon frequencies, magnon dispersions, static susceptibility and specific heat. 



PACS numbers: 75.10.Jm, 75.40. Cx, 75.40. Gb 
I. INTRODUCTION 

Vanadyl pyrophosphate (VOPO) presents a quantum 
magnetic system of S = | V 4+ ions, whose antiferromag- 
netic (AF) Heisenberg superexchange interactions giva 
a singlet ground state, with spin gap A = 3.1meV.EJ 
While early powder susceptibility and inelastic neutron 
scattering (INS) experiments were consistent with both 
spin-ladder and alternating-chain models.! 2 ] the first INS 
measurements on aligned single crystals^ confirmed the 
magnetic structure to be of alternating-chain conforma- 
tion. In conjunction with observations on the structurally 
related material VODPO4. iD20,Q the following picture 
was verified: the strongest exchange path (Ji) is the 
double V-O-P-O-V link through two phosphate groups 
along b; the next strongest ( J2) is the double V-O-V link 
between edge-sharing VO5 square pyramids, also along 
b; the structurally dimerized chains have coupling ratio 
A = J%j J\ ~ 0.8.U The V-O-V bond along a is very weak. 
These results are consistent with the single electron on 
V 4+ occuping the d xy orbital {be plane). 

The same experimenter^ obtained detailed measure- 
ments of a second, low-lying, triplet excitation with gap 
5.7meV and significant intensity over at least half of the 
Brillouin zone, while in addition the coupling between 
dimerized chains was found to be weakly ferromagnetic 
(FM). Interpretation of the latter features was offeredQ 
in terms of a frustrated coupling of the dimerized chains, 
via the long but presumably not insignificant V-O-P-O- 
P-O-V pathway, which would promote a triplet, two- 
magnon boundj-state. Although VOPO has a compli- 
cated structurejj which contains 8 V atoms per unit cell, 
the almost identical interatomic distances dyy suggested 
a model] treating all dimers as identical magnetic unit 
cells. 

Since the proposal of this picture, three important, 
additional experiments have been conducted on VOPO 
single crystals, namely further inelastic neutron scatter- 
ing (INS) measurements, Nuclear Magnetic Resonance 
(NMR) studies and Raman light scattering. Each has 



offered new, and sometimes surprising, additional infor- 
mation concerning the physics of VOPO. The purpose of 
this work is to review these results and provide a new, 
minimal model which encapsulates all of the observed ef- 
fects. We proceed in Sec. II with a discussion of INS and 
NMR results which lead to a two-plane description. In 
Sec. Ill we consider the strong magnetoelastic coupling 
observed by Raman scattering, and present a theoretical 
treatment of the coupled magnon-phonon system by the 
flow-equation method. The consequences of this coupling 
are illustrated for a variety of experimentally measure- 
able quantities in Sec. IV. Sec. V contains a summary 
and conclusions. 



II. NMR AND BILAYER STRUCTURE 

Contrary to the expectation of magnetically identical 
V 4+ ions, NMR measurements of the Knight Shift and 
relaxation rate 1/Tiu give clear evidence for two distinct 
species of 31 P and 51 V nuclei in VOPO. While results for 
both atomic species were not fully consistent, the pres- 
ence of separate magnetic environments is quite unam- 
biguous. An interpretation requires two distinct planes of 
coupled, dimerized chains, with different exchange con- 
stants (Ji, J%)a,Bi resulting in different spin gaps Aoa 
and A 0B . „ 

This deduction is consistent with the structure,!!! where 
8 inequivalent V atoms occupy 2 dimers per chain and 
2 chains per cell: despite the very small differences in 
V-V separations, these dimer units are clearly not iden- 
tical. The very weak c-axis coupling, observed in the 
INS dispersion, acts to isolate the differing planes of cou- 
pled chains. The two triplet branches measured by INS 
have the simple interpretation of two one-magnon dis- 
persion curves, one from each plane: this explains, why 
their intensities have very similar q-dependences.t3 The 
reported gaps A s = 35K and A 0j 4 = 68KH are in excel- 
lent agreement with the results of Ref. S. 

To answer the question of why the J values in the two 
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planes are so different, we note that the previous treat- 
ment of all dimers as identical magnetic unit cells was 
based on the fact that differences in dvv are 0(1O -3 ) 
of the unit cell dimensions. However, inspection of the 
structural parametersB reveals that differences in the lo- 
cations of P atoms, and thus in the interatomic P-V dis- 
tances dpv, are 0(1O~ 2 ). The PO4 group is the essen- 
tial element mediating the dominant superexchange path 
V-O-P-O-V. While quantitative accuracy is still lacking 
in ab mifio-superexchange calculations, particularly for 
V systems,£j the interaction magnitude is well known to 
have a very strong dependence on both interatomic spac- 
ings and bond angles. Thus one may expect a strong dif- 
ference between values J\a and Jib, but similar values 
of J 2 a and Jib ■ 
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FIG. 1. Schematic representation of the VOPO system; A 
and B denote the two decoupled types of plane. 

Retaining the basic framework of Ref. ^, we consider 
thus a model of two independent types of planes of al- 
ternating chains (Fig. [I]). We-i begin by fitting the two 
independent magnon branchesB to determine the appro- 
priate, effective superexchange constants. This is per- 
formed within a "static model" , by which is meant one 
with temperature-independent effective interactions, in 
contrast with the situation in Sec. IV where we will 
consider a model with phonon dynamics. rTjhe magnons 
are described by mobile triplet excitationsJi 2 ] and the fit- 
ting procedure differs from that presented!] for the single- 
plane model in two minor respects. First, as mentioned 
in that work, we fit the square of each dispersion curve, 
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because this smoother quantity gives a superior fit to the 
same order of expansion. The parameters iiy, for 7 = 
A, B, are determined to third order from the quantities 
Uj in Eq. (1) of Ref. g atO 

uqo = 1 + 7A 3 /32 + 2/£ + 3/i 2 _/i+/4 - 2Aa 
+2a 2 - (15A 2 a + 6Xa 2 - 12a 3 )/16 

uxo = -A - A 2 /2 + + 2a + 5A 3 /32 

-7A 2 a/16 - 9Aa 2 /8 + 3a 3 /4 

uoi = 2^t_ + 3fi 3 _ /4 - A 2 /i_ (2) 

u 2 Q = (A 3 - 2A 2 a + 4Aa 2 - 8a 3 )/16 



U02 = — pt_A*+/2 

Tin = — 3A 2 /i_/4 — a 2 /i_ + A/i_/i+/2 



Here A = J2/J1 (Fig. [[]), fx+ = {J a + Jb)/J\ and 
|U_ = (J — Jb)/Jx- The transverse part of the dis- 
persion w 7 (q) is very sensitive to the difference /i_ 7 of 
the interchain exchange constants, but rather insensitive 
to their sum /i+ 7 . The parameter a — Jf/Ji denotes a 
frustrating, next-neighbor interaction (Jy) along the spin 
chains. It is not present in our minimal model (Fig. |l|) 
and is not used in the fits presented in this section, but 
is included here for reference from Sec. IV, where such 
a term is generated within the dynamical model. The 
second difference from Ref. g is that we do not determine 
/i+ 7 from the condition on the Curie temperature de- 
duced from the magnetic susceptibility, as this procedure 
may no longer be applicable in the presence of phonons 
coupled strongly to the spin subsystem (below) . Instead 
we deduce the values /z+ 7 from [the maximum x, nax of 
the measured suseptibility curve. til 
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FIG. 2. Fitted dispersion relations for the two elementary 
magnons, (a) in chain direction at q a — and (b) transverse 
to chains at qt = 0. Data from Refs. ^ and [1^. 

The fit (Fig. ||) yields the following effective parame- 
ters: 

(Jia, J2A, J aA , J bA ) = (144.9K, 95.5K, 11.6K, 16.4K) 
(Jib, J 2 b, J aB , J bB ) = (122.9K,95.0K,13.4K,18.6K). (3) 

The dimerization parameters Xa = 0.659 and Ab = 0.772 
are rather different in each set of chains, but the values of 
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Ji are similar, as expected from the considerations above 
(V-O-V superexchange) . These static- model results for 
the chain couplings Jj, and J2 are very close to those of 
the NMR experiment Jj and are also in accord with values 
obtained recently_by a detailed fit based on uncoupled, 
dimerized chains£3 In comparison with the latter, the 
present procedure has the advantage of containing the 
true dimensionality of the magnon dispersion, so that no 
discrepancy occurs between a gap from a chain fit and a 
gap including interchain coupling. Furthermore, we find 
that non-crossing magnon dispersion curves have a better 
physical justification due to the similarity of J2 superex- 
change paths in both types of plane. Pressure-dependent 
measurements and comparative ab initio studies would 
be required to shed further light on the variation of su- 
perexchange parameters with atomic positions. 

In comparison with the single-plane model, the frus- 
trated interchain coupling is significantly weaker: in ob- 
vious notation, jj, a A = 0.080, /i a s = 0.109, = 0.113, 
HbB = 0.151. Any bound states of the elementary 
magnons can be expected to be at best only vegy weakly 
bound, and in addition to have very low weight.Q Another 
significant difference comes in calculating the static sus- 
ceptibility x{T) in the twOf-elane model: the fit to single- 
crystal [Susceptibility dataE-3 is much improved in terms 
of T max ld (see Fig. |J). This arises because of the contribu- 
tion of the higher-energy magnon band. Because the cal- 
culation of x(T) is approximate, in that it involves only 
one-magnon contributions with interaction effects on a 
mean-field level, remaining discrepancies between theory 
and experiment are not unexpected. However, this last 
point cannot account for the inconsistency observed in 
the high-temperature regime, where the moments should 
behave independently, and we return to this issue below. 

Two features remain the same as the previous frame- 
work. First, we contend that two, mutually frustrating 
AF interactions between the alternating chains remain 
the most likely scenario to account for the form of the 
FM a-axis dispersion. This is fitted appropriately by the 
weaker values (J a , Jb) above. Although this may only be 
treated qualitatively, the condition on the Curie temper- 
ature, 

B C w = -\ Yy ji -< + j2 -< + 2 ( ja -< + w 

7 

provides further support for frustrated couplings, because 
the chain couplings alone remain unable to satisfy this 
sum. Second, the two-dimensional dispersion relation 
in (q a ,Qb), with logarithmic singularities in the density 
of states, remains essential to explain the temperature 
scales of thermodynamic propertieSrSuch as xCO an d the 
electron spin resonance absorptionE 

III. SPIN-PHONON INTERACTIONS 

Recent Raman Jight-scattering experiments on single 
crystals of VOPOO have shown two important new fea- 



tures. One is a strong anharmonicity, in the form of 
a hardening (frequency increase) with decreasing tem- 
perature, of certain phonons, and the other a strong 
quasielastic scattering of magnetic origin. Both features 
are clear, qualitative evidence for magnetoelastic cou- 
pling. In the phonon system, the frequencies and polar- 
izations of strongly renormalizcd phonons indicate oscil- 
lations of the phosphate (PO4) groups. In particular, the 
123cm _1 phonon observed in aa polarization (in-plane 
motion transverse to the chains, expectedlij to have the 
strongest coupling) loses intensity very rapidly to the spin 
degrees of freedom. In the spin system, quasielastic scat- 
tering originates in energy fluctuations of the spins, which 
are strongly enhanced by temperature due to phonon 
coupling. In addition, the observed 2-magnon intensity 
falls rapidly with temperature; this feature seems to cor- 
respond to the onset of a continuum at u = 2Aqb, rather 
than to a singlet bound state as suggested in Ref. [IT]. 

At the qualitative level, this spin-phonon coupling is 
considerably stronger than that observed in the spin- 
Peierls system CuGeOs. We also note here that the 
authors of Ref. [17] do not observe magnetic scattering 
at u> = 2AoAn and state that this invalidates the 2- 
plane scenario.Q We suggest from i) the weakness of the 
w = 2Aos feature and ii) the low intensities associated 
with any calculation of bound states in the current pa- 
rameter regime, that the onset of the second continuum 
is simply too weak to have been detected here. The 
reported observation of this feature in further Raman 
experimental^ verifies this hypothesis. 

We wish to provide a theoretical description for the 
consequences of a significant magnetoelastic coupling, 
in terms of those phonons most strongly coupled to a 
spin system represented by the above, minimal model 
of triplet magnon excitations. Spin-phonon coupling ef- 
fects may readily be envisaged within a conventional, di- 
agrammatic treatment,L3 where perturbative inclusion of 
a magnon-phonon vertex with coefficient g would lead to 
0(g 2 ) self-energy corrections to the phonons due to the 
presence of magnons, and conversely. However, this type 
of approach would appear to be precluded here by the 
difficulties inherent in expressing a propagator for hard- 
core bosons, and in substituting frequency summations 
with constrained thermal occupation functions. This ap- 
proach has been applied-to consider spin-phonon coupling 
in a two-chain ladder,E2l by the introduction of an infi- 
nite but fictitious repulsion term for triplet excitations on 
the same bond. In addition to this weakness, the bond- 
operator method appears to be applicable at the quan- 
titative level only for very strongly dimerized systems. 
Here we will instead reproduce the mutual renormaliza- 
tion of phonons and magnons by a flow-equation method 
particularly suitable for local spin-phonon problems,Eil in 
which a unitary transformation is applied to the Hamil- 
tonian of the coupled system to eliminate the coupling 
term. 

Motivated by the qualitative observation concerning 
the importance of the PO4 groups, which is further sup- 
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ported by experimental observation on the phonon modes 
involved, in VOPO and related compounds ,E£I we will 
consider the previous model with phonon coupling only 
to the dominant dimer bonds J\a, Jib- We begin with 
the Hamiltonian in the form 



m=l,2 
m=l,2 

c 6t^ + G(6 J : + fe 4 t )S l 1 J ,S^}, 



+ 



where i and j are respectively indices for the dimer bonds 
( Ji) along and across the chains in each decoupled plane, 
and m = 1, 2 denotes the left or right spin in each dimer. 
The phonons {bi} are introduced as local, Einstein modes 
of fixed frequency loq at each dimer bond i; in reciprocal 
space these are nondispersive, and have the same weight 
at all wave vectors q. Because these phonon modes in- 
volve motion of the PO4 groups, and not of the mag- 
netic (V 4+ ) ions themselves, correlations between dis- 
placements may safely be neglected, and the approxima- 
tion of Einstein phonons justified. The spin-phonon cou- 
pling constant G is a free parameter to be fixed from ex- 
periment, but may in principle be very large: comparison 
with CuGeOs-suggests that values exceeding 0.3Jia are 
not excludedEj The term in Jf 1 describes a frustrating, 
next-neighbor coupling along the dimerized spin chains; 
this is zero in the bare model (Sec. II), but is generated 
at second order in G by the unitary transformation which 
eliminates the final term to leave only phonon terms bi- 
linear in {bi} in the resulting, effective Hamiltonian. A 
systematic discussion of the transformation procedure is 
presented in Ref. |2^. 

In deriving the effective Hamiltonian we retain only 
the leading order in G/luq, which is (G/wo) 2 ; omission 
of next-order terms (G/loq) 4 can be expected to be well 
justified. The transformation also involves an expansion 
in J/ojq, where terms 0(J/ujq) are retained but those 
0(( J/wo) 2 ) omitted; the validity of this approximation 
is not apparent in VOPO, which from the values of Ji 7 
and wo is not well in the anti-adiabatic limit, but is mo- 
tivated by the good, semi-quantitative agreement, with 
simulations for a similar model by Bursill et al.fa and 
may also be justified a posteriori. 

The general form of the effective Hamiltonian may be 
represented as 



H — H spin + H p h c 



AH X + AH Y + AH Z 



(6) 



in which H sp i n and i? p honon denote respectively the pure 
spin and phonon parts of Eq. (^). The coupling term is 
transformed into two correction terms in the spin sector, 
which following the notation of Ref. we write as 



and 



(7) 



(8) 



Here Aij den 
as C"S, _.S, 



ites the local coupling, which we have taken 
B and CA = [Hs,A] denotes the commu- 
tator of the quantity A with the spin-only part H$ of 
the starting Hamiltonian (Eq. (|B|)). A further correction 
term arises for the phonon sector, 



AH, 



(9) 



1,3 



where the expectation value (■ ■ -)s is computed for the 
spin sector. This last contribution was not considered 
in Ref. [23L where the focus was on the spin sector, and 
arises in the phonon sector from a mean-field treatment 
of the terms in Eq. (9e) of that work. 

Explicit evaluation of the additional terms generated 
by the transformation (Eqs. (R-0)) yields 



AH 



x 



r 2 

-— E( S .V S 



2 -) 2 



c 2 

km 



(10) 



to within a constant, using that for spins S = 1/2, 
(Si.S 2 +3/4) = Si.S 2 + 3/4 because the eigenvalues 
of the right-hand side are and 1. AHy and AHz are 
calculated from the result 



.1, ,. 



= \G 2 { J 2 (S i+1 j - 

— (Jo - Jb) (Sjj - 



s i-i,j) 



i) (11) 

(Sjj-i-Sij-i)] } > 



where four-spin terms are neglected. With due atten- 
tion to the number of neighboring chains, the additional 
terms AHx, AHy, and AHz may be transcribed into 
the form of the initial Hamiltonian (||), as corrections to 
the couplings {J}, 



AJi 7 = 
AJ 2l = 
AJ fl = 
AJ ai = 
AJ bl = 



2w 
G 2 



75- J2-Y c °th l — i 
2uji 



G 2 



„- J21 coth ( — -1 

4W| 



2lu 2 [ ai 



\2T ) 
\2TJ 



J b7 )coth^) 



2lo 2 {Jai 



Jb~/ ) coth 



(—) 

\2TJ 



(12a) 
(12b) 
(12c) 
(12d) 
(12e) 
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and as a correction to the phonon frequency in each plane 
of 

Ac 07 = %J 2l «S}.S} +1 ) - (Sf.S| +1 » . (13) 

The latter result was obtained by neglecting additional 
interchain coupling terms, on the grounds that both the 
coupling constants \J a — Jb\ ^ J% an d the spin ex- 
pectation values are much smaller between the chains 
than within them. The expectation values in Eq. (|l3| ) 
may be computed for independent chains by a high- 
temperature expansion in combination with the known 
T = behavior as described in App. A. 

These values are all finite at zero temperature. In order 
to perform a self-consistent calculation, we deduce the 
appropriate bare values of the coupling constants such 
that their renormalized values correspond to the exper- 
imental data. Thus the fitting parameters provided in 
the previous section are effective quantities already con- 
taining the low-T renormalization. We now discuss the 
quantitative properties of these corrections in the context 
of phonon anharmonicity, magnon thermal renormaliza- 
tion, and thermodynamic properties. 



with experiment appears realistic. For completeness, we 
add here that the results in Fig. || were computed with 
the bare coupling constants 

(Jia, J2A, JaA, JbA) = (132K, 112K, 10.9K, 17.1K) 
(Jib, J 2 b, Jas, J bB ) = (HIK, 110K, 12.5K, 19.5K). (14) 

Comparison with (g) shows that the coupling G has a 
strong effect on the chain superexchange parameters Ji 
and J2, but does not alter strongly the interchain values 
J a and Jb, whose sum we have determined as in Sec. II. 
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IV. EXPERIMENTAL CONSEQUENCES 



A. Phonon anharmonicity 



One observes from Eq. ( |13| ) that Aojq > 0, because 
for an AF system (Sf-S i+1 ) is negative. Further, one 
expects that Aluq — > as temperature becomes very large 
and the expectation values vanish. Both of these features 
are consistent with experiment. More interestingly still, 
Acjo does not vary monotonically with temperature, as 
(Sj.Sj) for nearest neighbors increases in magnitude at 
low temperatures; this result is exactly in accord with 
the small, initial increase observedEZI with temperature 
for both 70cm -1 and 123cm -1 phonons. 

Fig. H shows the thermal renormalization effect for 
an Einstein phonon of bare frequency loq — 118.4cm -1 
(14.68meV). The bare frequency value is chosen such 
that the effective frequency including corrections coin- 
cides with the experimental value of 123cm -1 at T = 0. 
The solid and dashed lines correspond to the different 
effects of the magnon modes in each type of plane for 
the same spin-phonon coupling. The coupling constant 
G = 6.6meV was chosen to reproduce the initial up- 
turn seen in cc polarization. Note that we do not wish 
to imply a connection between the two plane types and 
the two polarizations: each polarization measures some 
combination of the contributions from each plane type, 
determined by possibly different coupling constants to 
each. However, the results in Fig. [| provide a good in- 
dication of general magnctoelastic coupling effects, and 
of the splitting in energy of a single phonon mode when 
further couplings are neglected, which from comparison 



FIG. 3. Thermal renormalization due to magnons, for 
a phonon of frequency u)q = 118.4cm -1 with coupling 
G — 6.6meV. Symbols are experimental data in cc and bb 
polarizations, from Ref. [l?]. 

The most obvious qualitative feature in Fig. ^ is that 
magnon effects alone cannot account for the observed 
anharmonicity over the full temperature range. This ap- 
pears to consist of two components: at low temperatures, 
magnon-related renormalization effects occur on the en- 
ergy scale of the magnon gap Aqb — 35K; at higher tem- 
peratures, a further decrease of ojq is present, on the 
energy scale (180K) of ujq itself. The latter behavior 
we ascribe to the lattice structure, and await dynami- 
cal simulations, measurements of the thermal expansion 
coefficients, and magnetostriction-experiments to verify 
this point. Recent measurementsEa of ultrasound prop- 
agation in (VO)2P207 performed in high, pulsed fields 
found a strong magnetoacoustic coupling, and indicate 
that the required lattice information may be readily ex- 
tracted. The small maximum in the phonon frequency at 
low T would not be expected from lattice anharmonic- 
ity, and the explanation in terms of a magnon renormal- 
ization justifies its use in fitting the coupling constant. 
The resulting value of 0.58 for g = G/J\a {Jia taken 
from Eq. (14)) reflects the very strong coupling present 
in VOPO. It is not implausible in view of the complicated 
superexchange path and its demonstrated sensitivity to 
the atomic position of P, and not inconsistent with expec- 
tations based on a comparison with CuGeC>3. We note 
finally that all q-dependence of the phonon modes has 
been omitted from our discussion, for the reasons speci- 
fied in the previous section. 
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B. Magnon renormalization 

Inspection of the renormalized coupling constants 
(Eqs. (|l2[)) reveals the following qualitative features, i) 
the strong interaction constant J\ has a temperature- 
independent enhancement; ii) the weaker bond J2 is sup- 
pressed with increasing temperature; Hi) in-chain frus- 
tration Jf is introduced, and rises with increasing T; iv) 
the interchain coupling diminishes as temperature rises. 
Based on these observations, one expects a band flatten- 
ing in the chain direction: the average frequency depends 
mostly on Ji, and is little changed, whereas the differ- 
ence between the extrema, w max — w m ; n depends mostly 
on J2 (and partly on Jf), and decreases with increasing 
temperature. 



diction: INS measurements could be expected to reveal 
large upward shifts in the magnon dispersion peak po- 
sition at the zone centre, and downward shifts at the 
zone boundaries, even at temperatures where broaden- 
ing of the signal does not compromise the fit quality. 
Our prediction for this effect is shown as a function of 
temperature in Fig. ^(b). Finally, we note that similar, 
very strong band-flattening features have been observed 
in INS measurements performed over a range of tem- 
peratures for the material KCUCI3, where complicated 
superexchange paths are also known— to be important in 
determining the magnetic couplingsH3 

C. Thermodynamic quantities 
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FIG. 4. (a) Thermal (phonon) renormalization of ele- 
mentary magnon modes at q a = 0, illustrated for an Ein- 
stein phonon with frequency luq — 118.4cm -1 and cou- 
pling G = 6.6meV, at temperatures of 10K and 300K. (b) 
Band-flattening effect of thermal renormalization, shown by 
thermal evolution of Auj — lu(0, n) — u(0,0), the difference 
between dispersion extrema. 

Fig. ||(a) illustrates this effect by comparing the 
magnon dispersion relations in the chain direction, de- 
termined at low and high temperatures for a spin system 
coupled to an Einstein phonon mode at luq = 118.4cm -1 , 
with the same spin-phonon coupling, G — 6.6meV, deter- 
mined from the anharmonicity. Flattening of the magnon 
bands as a consequence of this phonon coupling is clearly 
very strong. This result leads to a readily falsifiable pre- 



The magnon thermal renormalization arising as a re- 
sult of spin-phonon coupling should be manifest in cor- 
rections to thermodynamic quantities at high tempera^, 
tures. Considering first the static susceptibility x(T),L3 
this may be computedQ by integration over the available 
magnon modes. x(T) grows initially with temperature, 
its behavior governed exponentially by the gaps, has a 
maximum on the order of the dominant energy scales J\ , 
and converges at high-T to Curie- Weiss behavior. Fig. ^ 
compares the single- crystal experimental susceptibility 
with the single-plane model predictionsQ, the two-plane, 
static model for the parameters fitted in Eqs. (||), and 
the two-plane, dynamical model with thermally renor- 
malized magnon dispersions (Eqs. (pr|)). In the last case, 
one Einstein phonon mode per bond (ujq = 118.4cm -1 , 
G = 6.6meV) is used, as in the calculations of the previ- 
ous section. 
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FIG. 5. Comparison of experimental susceptibility x{T) 
(circles) with single-plane model (dashed), static, two-plane 
model (dot-dashed), and dynamical, two-plane model con- 
taining renormalizing effects of an Einstein phonon with fre- 
quency uiq — 118.4cm -1 and coupling G = 6.6meV (solid). 
Experimental data from Ref. 
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As noted in Sec. II, the two-plane, static model gives a 
rather good account of the low-temperature regime and 
the maximum, and remains quite satisfactory at higher 
temperatures. It is perhaps surprising in view of the 
large coupling constant G that the effects of dynamical 
phonons are so small. However, it is also evident that 
the model including thermal renormalization of the su- 
perexchange parameters provides excellent agreement in 
all three temperature regimes, including high-T. This re- 
sult not only confirms the consistency of the method, and 
of the value of G deduced from the phonon anharmonic- 
ity, but has a profound consequence for interpretation 
of high-T susceptibility data in terms of a Curie- Weiss 
temperature in the presence of phonons. The effects of 
a spin-phonon coupling on Gcw may be computed rig- 
orously, as shown in App. B. It is found that the value 
of 0cw including spin-phonon coupling corresponds to 
that given in Eq. (Q) if the renormalization of Ji 7 due 
to AHx (0) is taken into account. Corrections due to 
AHy @ are found to cancel. We note finally that in the 
real material, several types of phonon mode contribute to 
a susceptibility renormalization. However, because their 
effects are proportional to the square of the correspond- 
ing coupling constant, contributions beyond that from 
the dominant phonon mode may be weak. 
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FIG. 6. Specific heat of magnetic system, compar- 
ing static (solid) and dynamical (dashed) two-plane models. 
Dot-dashed line is for dynamical model in a magnetic field 
of 24T. Points are experimental data from Ref. ^| Dy- 
namical models contain one Einstein phonon of frequency 
u>o = 118.4cm -1 and coupling G = 6.6meV. 

We close this section by computing a further exper- 
imentally accessible quantity which is expected to illus- 
trate thermal renormalization due to magnetoelastic cou- 
pling. The specific heat of the spin system alone is diffi- 
cult to isolate from phonon contributions, and to date has 
been deduced only indirectly from quasielastic Raman 
scattering. EZl The calculation proceeds as for the mag- 
netic susceptibilityQ but with integration in reciprocal 



space over the magnon contributions to the second tem- 
perature derivative of the free energy, and may be readily 
extended to finite magnetic fields. The results presented 
in Fig. [| demonstrate again both the utility of the static 
model and the small but significant upward renormaliza- 
tion offered by the dynamical model at high T. -Both 
sets of data are close to the results of calculationsE3 for 
a single dimerized chain with appropriate A, indicating 
that interchain coupling effects have only a small effect 
on the overall features of thermodynamic quantities. By 
contrast, a strong magnetic field is most effective at low 
T, where it acts to reduce the gap, and at intermediate T 
where the peak value of C(T) is suppressed. The exper- 
imental data show too much scatter to be of particular 
utility, and we await more refined measurements of this 
quantity. 

V. SUMMARY 

In conclusion, we have analyzed a two-plane model for 
VOPO, which is fully consistent with the known crystal 
structure, and emphasize the key role of the PO4 groups 
mediating the dominant superexchange path. A full un- 
derstanding of the magnetic properties is not possible 
without considering the strong magnetoelastic coupling 
effects observed in Raman light-scattering experiments. 
We account for these qualitatively in terms of phonon 
modes of the phosphate groups, and quantitatively by 
using the flow-equation method. Phonon hardening is 
found to have contributions from both spin coupling and 
lattice anharmonicity, although only the former can ac- 
count for the observed low-temperature softening, which 
provides an estimate of the spin-phonon coupling con- 
stants. 

We obtain an excellent fit to the magnon dispersion 
data for the two triplet modes observed by INS, and 
find complete agreement with the parameters deduced 
from NMR studies. In addition, there is a strong ther- 
mal renormalization of the magnon dispersions due to 
phonons, in the form of a band flattening, which should 
be clearly visible in INS measurements performed as a 
function of temperature. We show further that, while 
the static two-plane model parameters give a good ac- 
count of static susceptibility and specific heat data, the 
renormalization effects of dynamical phonons lead to ad- 
ditional corrections which are required to reproduce the 
high-temperature limit. However, we have found that 
even for very strong coupling, the influence of phonons 
on thermodynamic magnetic properties is rather small, 
and such quantities may be quite well described by a 
static model. Unambiguous evidence of phonon effects is 
provided only by dynamical properties, such as Raman 
spectra and triplet dispersions. 
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APPENDIX A: DLOG-PADE EXTRAPOLATION 
OF EXPECTATION VALUES 

In this appendix we present details of the computa- 
tion of expectation values such as (Sj --Sf j), as func- 
tions of temperature. The basic physical consideration 
is straightforward for a gapped spin system, and we will 
demonstrate it for the expectation value A(T) = (Sj-Sf) 
of the spin correlation function for the strong bond in 
a dimerized spin chain (i is the dimer index). At zero 
temperature, A(T) has the finite, negative value Aq. At 
small but finite temperature, A{T) will deviate from Aq, 
with an exponentially small deviation due to the presence 
of the spin gap. It is plausible to take this deviation to be 
positive, A(T) — Aq > 0, and so one expects (c/. Fig. [?]) 

A(T) w A + A\T V exp(-A/T) for T small (Ala) 
A(T) oc 1/T for T large, (Alb) 

where A is the spin gap. The exponent v = d/2 is deter- 
mined for a dispersion with quadratic minima only by the 
dimensionality d. For the quantitative computation we 
use two sources of input, namely i) an expansion at zero 
temperature about the dimer limit, such as that given in 
Ref. and ii) a high temperature expansion such as 
that presented in Refs. ^6 27. 



From (i) one may deduce the values of Aq and A. For 
a sufficiently dimerized system such as (VO)2P20j, the 
plain truncated series can be used. The value of Aq is ob- 
tained from the expansion of the ground state energy as 
a function of the strong-bond coupling J\ , the weak-bond 
coupling J2, and the next-nearest-neighbor coupling Jf, 
by partial differentiation with respect to J\ . 

The extrapolation is performed using a Dlog-Pade ap- 
proximant for A(T) — Aq in the following manner. Ex- 
pansion in (3 up to 9 th order is used to describe 



a:=d A(l/P)/(A(l//3)-A o ). 



(A2) 



Then, because we expect regular behavior in the limit 
P — > co, it is appropriate to change to variable u = 
13/(1 + (3) (whence (3 — u/(l — u)). Because u oc (3 for 
small values of u and /3, an expansion of a in terms of 
u is available up to order 8. This is used to set up a 
[8, 2] Pade approximant Ps/Q2 in u, where the subscripts 
stand for the degree of the polynomial in u. The miss- 
ing constants for order 9 and 10 are determined from the 
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known asymptotic behavior. From Eq. (Ala) we have 
a = -A - d/{20) + £>(/r 2 ), which yields 



APPENDIX B: CURIE- WEISS TEMPERATURE 
AND SPIN-PHONON COUPLING 



5 U (P 8 /Q 2 ) = 1/2 

u—1 



(A3) 
(A4) 



for chain geometry. In this way, reliable and non- 
defective interpolations between the zero-temperature re- 
sult and the high-temperature behavior can be obtained, 
as illustrated in Fig. 



i 




FIG. 7. Expectation value of strong-bond spin correlation 
function (Sl-S|), determined as function of temperature for 
J2/J1 = 0.75 and Jf/Ji = 0.05 by procedure described in 
text. 

Applying the above procedure to the phonon frequency 
shift ([[3]), one encounters a robust pole in the Pade ap- 
proximants. This pole is approximately independent of 
the actual Pade approximant, so is not spurious. In fact 
the procedure explained above must fail if A(T) — A 
changes sign, an occurrence which may not be excluded 
for general expectation values. Indeed, we find that the 
weak-bond expectation value in a dimerized chain does 
display such a sign change, although the difficulties this 
introduces are readily avoided. The function considered 
is separated into two parts, each without sign change. In 
practice, this was effected by adding to the expectation 
values in Eq. ( |l3| ) four times the strong-bond expectation 
value, and applying the above procedure to the combined 
expression. The strong-bond contribution from separate 
computation (above) is then subtracted to yield the final 
result. 

This technique, which relies on a physically appropri- 
ate combination of information from zero and high tem- 
peratures, permits rapid computation of the frequency 
shift despite the need to determine the effective couplings 
anew at each temperature. 



In this appendix we derive the effect of a magnetoe- 
lastic coupling on the Curie- Weiss temperature, Ocw- 
Because at high temperatures one assumes x(T) oc (T — 
©cw) _1 , this procedure involves computing the sublead- 
ing 1/T behavior in (T - Ocw) -1 = 1/T + 6 C w/T 2 + 
0(T~ 3 ). We consider first a uniform spin chain coupled 
to Einstein phonons, 



# = J ]T SrS l+ i(l + g{h + &})) + wb% 



(Bl) 



Direct expansion in 1/T of the partition function 
Tr[exp(— (3H)] is not possible due to the infinite- 
dimensional bosonic Hilbert space, and we employ in- 
stead an expansion in J. 
The first two terms read 

Tr[exp(-0H)] = 0( J 3 ) + ■ 
(l + 7VJ 2 Tr spinsl , 2 [(S r S 2 ) 2 ](/3 2 /2! + .g 2 /i/^o)) , (B2) 

with Zq = (1 — e _/3 ") _1 the bosonic partition function 
and the coefficient /1 given by 

,/3 ,/3 r /3 



fl 



JO 



e -/3 lW 6 T 6 Me -/3 2W b r fc ue -/33^b T 6 



, (B3) 



where u denotes b + tf. The coefficients f m are the gener- 
alization of (B3) to 2m /3-ordered factors of u. The eval- 
uation of Tr sp i ns i : 2(Si-S2) 2 is straightforward, and yields 
3/16. The evaluation of / is possible by explicit calcula- 
tion, or more easily by the observation that / equals the 
second coefficient in A in an expansion of Z(X), where 



Z(X) = z + J2 a 2 "7„ 



(B4) 



is the partition function of H = cub^b + \u. Transform- 
ing to shifted bosons b = b + X/ui, the shifted parti- 
tion sum is easily found to be Z(X) — Z exp(/3A 2 /w), 
whence it follows directly that f m = Zo(/3/uj) m /m\. Be- 
cause f m is of order (3 m Zo, the general phononic con- 
tribution in Eq. (B2) is of order J 2m ((3/uj) m , or higher 
in p. Hence one obtains a systematic high-temperature 
expansion circumventing the problem of the infinite di- 
mensional bosonic Hilbert space, which is cut off by the 



factor e 



-/3ub f b 



In analogy to the expansion (|B2|), the variance of the 
magnetization M := Sf is given by 



lTr[M 2 exp(-/3i0] = ifli 



R 2 + 0(J 3 ), (B5) 



where R\ is the partition sum (B2) and 







R 2 = 2Z^ (-/3JTr spinsl , 2 [S^(S r S 2 )] + 

J 2 Tr spinsl , 2 [^^(S r S 2 ) 2 ] (l3 2 /2 + g 2 f 1 /Z )) . (B6) 



The susceptibility x(T) is the ratio of Eqs. ( |B5| ) and ( JB2| ) 
to order J 2 , 

4T X = 1 - ^ - Ff. - + 0(J 3 f3 2 ) (B7) 

= l-/?(j + ^)+0(/3 2 ), (B8) 
from which follows 

This result implies quite generally that the Curie- Weiss 
temperature is lowered by spin-phonon interactions, i.e. 
in AF systems the modulus of &cw increases. For Ein- 
stein phonons the contribution of spin-phonon coupling 
to the Curie- Weiss temperature amounts to — G 2 /(8u>) 
per bond linked to each site with coupling G = gJ. This 
is exactly the same effect as the renormalization con- 
tained in AHx (Eq. ((?])), a fact which corroborates the 
validity and utility of the flow-equation approach. 

To conclude this appendix, we have for the model (j^) 
a change of Ocw due to spin-phonon interactions given 
by 

©CW;total = @CW;barc — G 2 /(8w) (BIO) 

where 0cw;bare is the value given in Eq. (^). The expan- 
sions of Eqs. (B2,B5|) provide an interesting and system- 



atic extension to higher orders in J. Generally, a calcu- 
lation in J 2m will provide results for a high-temperature 
expansion up to j3 m . 
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